Gamma-ray Activity in the Crab Nebula: 
The Exceptional Flare of April 2011 
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ABSTRACT 

The Large Area Telescope on board the Fermi satellite observed a gamma-ray flare in the 
Crab nebula lasting for approximately nine days in April of 2011. The source, which at optical 
wavelengths has a size of «11 ly across, doubled its gamma-ray flux within eight hours. The 
peak photon flux was (186 ±6) x 10~ 7 cm -2 s _1 above 100 MeV, which corresponds to a 30-fold 
increase compared to the average value. During the flare, a new component emerged in the 
spectral energy distribution, which peaked at an energy of (375 ± 26) MeV at flare maximum. 
The observations imply that the emission region was relativistically beamed toward us and that 
variations in its motion are responsible for the observed spectral variability. 
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1. Introduction 

The Crab nebula is the remnant of a supernova 
observed in 1054 AD. The explosion left behind 
a rotating neutron star emitting electromagnetic 
radiation pulsed at the rotation period, that pow- 
ers a wind of relativistic particles. These particles 
interact with the remnant gas and magnetic field, 
causing the nebula to glow brightly at all wave- 
lengths, predominantly by synchrotron radiation. 
The spectral energy distribution (SED) of the neb- 
ula is, accordingly, dominated by a synchrotron 
component, extending from radio wavelengths into 
the gainma-rav band (Hester 2008). Above 450 
MeV a second component emerges, attributed to 
inverse-Compton scattering by the same relativis- 
tic particles (Gould 1965; de Jager & Harding 
1992: Atoyan et al. 1996). The angular size of 
the Crab nebula is «0.1° in the optical and smaller 
at higher energies. This corresponds to 3.5 pc, or 
11 ly, at its estimated distance of 2 kpc (Trimble 
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1973). 

Today, the pulsar and nebula (henceforth re- 
ferred to together as the Crab) are considered 
prime examples of non-thermal sources in the Uni- 
verse and serve as a laboratory for relativistic 
plasma physics. New puzzles for our understand- 
ing of the Crab have been posed by the detection 
of three bright gamma-ray flares by the AGILE 
satellite and the Large Area Telescope (LAT) on 
board the Fermi satellite between 2007 and 2010 
(Abdo et al. 2011a; Tavani et al. 2011). Dur- 
ing these flares the unpulsed component of the 
gamma-ray flux increased by a factor of wlO on 
time scales as short as 12 hours (Balbo et al. 
2011), while the period and flux of the pulsed com- 
ponent remained stable. 

More recently, in April of 2011 the LAT de- 
tected a fourth flare, three times brighter than 
any of the previous ones (Buehler et al. 2011; 
Hays et al. 2011). The flare was swiftly con- 
firmed by the AGILE satellite (Striani et al. 
2011). Observations at lower energies, in partic- 
ular by the Chandra X-ray observatory, have not 
yet revealed any variability correlated with the 
gamma-ray flares (Tennant et al. 2011). How- 
ever, analysis of these observations is ongoing and 
will be discussed elsewhere. Here we present the 
LAT gamma-ray results obtained during the flare, 
together with a broader analysis of the first 35 
months of Crab observations by Fermi. 

2. The Large Area Telescope and data 
analysis 

The LAT is a pair-conversion telescope, sensi- 
tive to gamma rays with energies greater than 20 
MeV. It has a large field of view («2.4 sr) and 
images the full sky every three hours. The angu- 
lar resolution of the LAT varies with photon en- 
ergy. The 68% containment radius ranges from 
approximately 6° at 70 MeV to 0.2° above 10 
GeV (Atwood et al. 2009). LAT scientific obser- 
vations began in August 2008. We analyzed data 
taken within 20° of the Crab in the first 35 months 
of observations (MJD 54683-55728). The average 
spectral properties of the Crab were derived from 
the first 33 months of observations (MJD 54683- 
55664), excluding the April 2011 flare. 

Fluxes and spectra were obtained by maxi- 
mizing the likelihood of source models using un- 


binned gtlike from the Fermi Science Tools 9- 
23-01. The models included all sources in the 
second LAT source catalog within 20° of the 
Crab position (Abdo et al. 2011c) plus models 
for the Galactic and isotropic diffuse emission 
(gal_2yearp7v6jv0 , iso_p7v6source). The pa- 
rameters left free to vary in the likelihood fit 
for 33-month average spectra were the spectral 
parameters of the Crab, the normalization of 
the diffuse components and a power-law spec- 
tral index for scaling the Galactic diffuse emis- 
sion model. For analysis on shorter time scales 
only the isotropic diffuse normalization was var- 
ied along with the Crab spectral properties. All 
other parameters were fixed to the 33-month av- 
erage maximum likelihood values. The Sun was 
included in the source model when it was within 
20° of the Crab. The solar spectra during the 
two passages in front of the Crab were taken from 
Abdo et al. (2011b). The Moon w'as not included 
in the source model as its gamma-ray contamina- 
tion was found to be negligible for observations 
presented here. 

We used the P7_V6_SOURCE instrument re- 
sponse functions without in-flight PSF corrections, 
selecting photon events between 70 MeV and 300 
GeV. Compared to the more typical 100 MeV 
threshold this choice leads to additional system- 
atic errors due to increased dispersion in the pho- 
ton energy reconstruction. The overall systematic 
flux error is energy dependent: it amounts to 30% 
at 70 MeV and decreases to 10% above 10 GeV 
(Rando et al. 2011). The dominant part of this 
systematic error is related to the overall flux nor- 
malization. It is caused by uncertainties in the 
effective area determination and of the overall nor- 
malization of energy scale. 

The stability of the LAT instrument over time 
w'as tested for all time scales addressed in this 
publication using the Vela and Geminga pulsars, 
which are found to be stable in flux (Abdo et al. 
2010b, c). Their flux variations were <10%, yield- 
ing an upper limit on the variations of the sys- 
tematic errors with time. A more detailed study 
of the systematic uncertainties is currently' being 
prepared within the LAT collaboration for publi- 
cation. 

The pulsar phase was assigned to the detected 
gamma rays based on a high-time-resolution pul- 
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Fig. 1. — Phase profile of gamma rays above 70 
MeV within 3° of the Crab Pulsar for the first 33 
months of Fermi observations (black histogram) 
and during the April 2011 gamma- ray flare be- 
tween MJD 55663.70-55671.02 (black markers 
with error bars). The gray region indicates the 
adopted off-pulse interval where the emission is 
dominated by the nebula. The flare phase profile 
has been multiplied by a factor 59, such that the 
excess above the off-pulse counts is the same as for 
the 33 months of observations. This demonstrates 
that the flare is a phase-independent flux increase. 


sar ephemeris 1 . To obtain the latter we extracted 
400 pulsar times-of-arrival (TOAs) from LAT pho- 
tons collected from MJD 54684-55668 with a typi- 
cal uncertainty of ~ lOOps (Ray et al. 2011). Us- 
ing tempo2 (Hobbs et al. 2006), we fit a timing 
solution to these TOAs with a typical residual of 
108/iS, or about 3 x 10~ 3 of the pulsar period. To 
obtain these white residuals, we modeled the pul- 
sar timing noise using the method of Hobbs et al. 
(2004), with 20 harmonically related sinusoidal 
terms. 

3. Time-Average Energy Spectra 

The Crab appears to the LAT as a point source, 
even at the highest photon energies. To sepa- 
rate nebular gamma-ray emission from that of the 
pulsar we apply a cut on the pulsar phase. The 
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Fig. 2. — Spectral energy distribution for the 
Crab nebula averaged over the first 33 months of 
Fermi observations. The axis on the right side 
indicates the isotropic luminosity. Also shown 
are data from COMPTEL in the soft gamma- 
ray band (Kuiper et al. 2001) and very high en- 
ergy gamma-ray measurements from Cherenkov 
telescopes (Aharonian et al. 2006; Albert et al. 
2011). The dashed line shows the maximum like- 
lihood model in the parametrization described in 
the text. 


phased count rate of the pulsar with its double- 
peaked structure is shown in Figure 1. The pul- 
sar dominates the phase-averaged gamma-ray flux, 
but its flux in the off-pulse interval from 0.535 to 
0.885 is negligible (Abdo et al. 2010a). It is in 
this interval that we measure the properties of the 
nebula. 

The LAT detects the nebula in the energy 
range between the the high-energy end of the syn- 
chrotron and the low-energy end of the inverse- 
Compton components of the SED. The average 
nebular spectrum measured during the first 33 
months of Fermi observations is shown in Fig- 
ure 2. We fitted it as the sum of synchrotron 
and inverse- Compton components. The differen- 
tial photon spectrum, $(E), of the synchrotron 
component was parametrized with a power law 


<M£) = 


F s (is - 1) 
(100 MeV) 1 ' 


E~ ls , 


( 1 ) 


A The ephemeris parameter file will be made available at 
http: / / fermi . gsfc .nasa. gov /ssc/ data/ access /lat/ ephems/ 
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with an integral photon flux above 100 MeV of 
Fs = (6.1±0.2) x 10~ 7 cm“ 2 s“ x and a soft photon 
index of -ys — 3.59 ± 0.07. 





The spectrum of the inverse-Compton compo- 
nent softens significantly with respect to a power- 
law at higher energies, as expected from mea- 
surements at very high energy (VHE) gamma 
rays (Aharonian et al. 2006; Albert et al. 2011; 
Abdo, et al. 2010a). The spectrum was parametrized 
by a smoothly broken power law with a curvature 
index of p = 0.2: 


$i(E) = <hj, 0 



( 2 ) 


The spectrum has a photon index of 7/4 =1.48 
± 0.07 at low energies and softens to 7^2 = 
2.19 ± 0.17 above a break energy of Eb = (13.9 
± 5.8) GeV. The flux normalization is ‘L/o = 
(5.6 ± 1.4) x 10 -10 cm -2 s -1 MeV -1 and the inte- 
gral flux above 100 MeV is F / = (1.1 ±0.1) x 10 - ' 
cm -2 s -1 . 

The averaged pulsar spectrum in the first 33 
months of observations was measured in the on- 
pulse after accounting for the nebula emission. We 
parametrized the pulsar spectrum with a power- 
law function with a super-exponential cutoff 

*P(E) = Dp.o Vter (3) 

The best-fit value for the spectral index is 7 p = 
i.59±0.01, for the break energy Ep tC = (504±63) 
MeV and for the curvature index k — 0.43 ± 0.01. 
The normalization is given by $p, 0 = (8.1 ±0.5) x 
10 -10 cm -2 s -1 MeV -1 . The integral flux above 
100 MeV is Fp = (20.4 ± 0.1) x 10 -7 cm -2 s -1 , 
in agreement with the value measured after eight, 
months of observations (Abdo et al. 2010a). 

The pulsed emission from the Crab has recently 
been detected between «25-400 GeV. This emis- 
sion is constrained to pulses which are narrower 
by a factor of approximately two compared to the 
LAT energy range. Pulsed emission at these ener- 
gies challenges current pulsar models (Aliu et al. 
2008. 2011; Aleksic et al. 2011). A more detailed 
spect ral analysis of the Crab pulsar in the Fermi- 
LAT energy range and its connection to the pulsed 
VHE emission is currently being performed and 
will be published elsewhere. 


4. Temporal flux variations 

In order to probe the flux variation over time, 
the flux from the Crab was evaluated in 12-hour 
time intervals. The combined pulsar and nebular 
spectrum was modeled as a power-law in energy, 
as further spectral features are not resolvable on 
these short time scales for typical observed fluxes. 

The resulting light curve is shown in Figure 3, 
where the three flares previously reported by the 
Fermi-LAT collaboration are indicated. Variabil- 
ity is seen over the whole 35-month period. This 
is also apparent in the Fourier power density spec- 
trum (PDS) shown in Figure 4, which is approx- 
imately described by a power law with a spectral 
index of «0.9. The PDS was derived using the fast 
Fourier transform for evenly spaced data, inter- 
polating over the short time interval MJD 54901- 
54906, during which the Crab was not observed by 
Fermi We verified that results are essentially in- 
distinguishable from power spectra computed us- 
ing the techniques for unevenly sampled time se- 
ries described later in this section. The PDS shows 
significant power above the noise level on time 
scales from years to weeks. Variations on these 
time scales are also present outside of the flaring 
periods. This can be seen in the PDS of the time 
interval MJD 54884-55457 also shown in Figure 4, 
where the nebula did not show large variations in 
flux. The spectral index of the PDS during this 
time is approximately 1.0. 

There are compelling reasons to believe that 
the phase-averaged Crab pulsar and the inverse- 
Compton component flux are constant in the 
LAT energy range (de Jager & Harding 1992; 
Atoyan et al. 1996), or at least that they vary 
only very slowly as the pulsar spins down, and 
that any rapid variability is attributable to the 
synchrotron component of the nebula. This was 
confirmed observationally on monthly time scales 
for the first 24 months of Fermi observations 
(Abdo et al. 2011a). To test this assumption on 
the 12- hour time scales studied here, we measured 
the flux in the on-pulse interval, subtracting the 
nebula flux measured in the off-pulse interval for 
each time bin. The pulsar flux was found to be 
stable within 20%. No significant variations were 
found. The inverse-Compton component of the 
nebula is too faint to be consistently detected in 
12-hour time windows. Therefore the variability 
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Fig. 3. — Integral flux above 100 MeV from the Crab as a function of time for the first 35 months of Fermi 
observations. The time binning is 12 hours. The dotted blue line indicates the sum of the 33-month average 
fluxes from the inverse-Compton nebula and the pulsar. The dashed blue line shows the average flux of 
the synchrotron nebula summed to the latter. Fluxes are shown on a log scale. The three flares detected 
by the Fermi LAT are indicated. A fourth flare was detected in October 2007 (MJD 54380) by AGILE 
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Fig. 4. — Crab nebula Fourier power density spec- 
trum (PDS), calculated from the light curve of 
the first 35 months of Fermi observations shown 
in Figure 3. The PDS of the full time interval 
is shown by the solid green line (scaled down by 
1/100 for better visibility). The PDS of the low ac- 
tivity period between MJD 54884 55457 is shown 
by the solid blue line. The PDS of the April 2011 
flare is indicated by the solid red line and was 
calculated from the light curve shown in Figure 
5. A smoothing with a running average of four 
bins was applied to all spectra. The PDSs ob- 
tained before smoothing are shown in colored dot- 
ted lines. Black lines show the best fit function of a 
power-law function (dashed) plus a constant white 
noise component (solid) for the unsmoothed spec- 
tra. The best fit spectral indices are given in the 
text. Dotted black lines indicate the ±ler, the +2cr 
and +3 a confidence intervals derived from white 
noise simulations for the April 2011 flare PDS. 

of this component cannot be quantified from our 
analysis on these time scales. However, significant 
flux variations of the inverse-Compton component 
on short time scales are strongly disfavored the- 
oretically. The observed absence of variations on 
monthly time scales is consistent with constancy 
on the shorter scales of interest here. 

On 2011 April 9 the flux of the Crab increased 
dramatically. Once the source reached flux levels 


comparable to the 2009 and 2010 flares the Fermi 
satellite was commanded to switch to a pointed 
observation targeting the Crab (MJD 55663.70- 
55671.02). In this mode the exposure toward the 
source increased by a factor of about four com- 
pared to the standard all-sky monitoring. During 
these observations the Crab erupted to a peak flux 
of (186 ±6) x 10~ 7 cm -2 s _1 above 100 MeV dur- 
ing the 12-hour period centered at MJD 55667.14, 
as shown Figure 3. This corresponds to a flux 
increase of a factor 7 compared to the average to- 
tal flux from the Crab and a factor of 30 com- 
pared to the flux from the nebula. We verified that 
the flare is indeed positionally coincident with the 
Crab nebula, with a best-fit localization of E. A. = 
83.65° Dec. = 21.98° (J2000), and an error radius 
of 0.04° at 95% confidence. 

The high gamma-ray flux during the flare, com- 
bined with the increased exposure of the pointed 
observations allowed us to study the flux evolution 
down to time scales below the wl.5 hour orbital 
period of Fermi. On such time scales occupation 
of the source bv the Earth needs to be taken into 
account. The time binning was therefore adjusted 
such that only time periods during which the Crab 
is visible to the LAT are used. These visibility win- 
dows were further split into bins of equal exposure, 
yielding a mean bin duration of nine minutes. The 
evolution of the flux during the flare in this binning 
is shown in Figure 5. The flare lasted for approx- 
imately 9 days and is composed of two sub-flares, 
peaking around MJD 55665 and MJD 55667. Dur- 
ing both, the flux increased rapidly, reaching its 
maximum value within approximately one day. A 
second rise or “shoulder” is observed during the 
decaying phase of both sub-flares. Whether this 
is coincidental cannot be assessed on the basis of 
these two events alone. If these shoulders are not 
interpreted as additional flares of lower amplitude, 
the decay time is approximately 1.5 and 3 days for 
the first and second sub-flare, respectively. 

To investigate the significance of peaks on 
smaller time scales in the flare light curve, we 
decomposed it into periods compatible with a con- 
stant flux. The time windows are referred to as 
“Bayesian Blocks” (BBs). They were determined 
finding the partition which maximizes the sum of 
the cost function assigned to each BB (Seargle 
1998). The cost function was set to the loga- 
rithm of the maximum likelihood for the constant 
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Fig. 5. — Integral flux above 100 MeV as a function of time during the 2011 April Crab flare. The light 
curve is binned into equal exposure bins during times with no Earth occultation, with a mean bin duration 
of nine minutes. The dotted line indicates the sum of the 33-month average fluxes from the inverse- Compton 
nebula and the pulsar. The dashed line shows the flux of the average synchrotron nebula summed to the 
latter. The solid black lines show the best fit of a model consisting of a constant plus an exponential function 
at the rise of both sub-flares (see text). The blue vertical lines indicate the intervals of each Bayesian Block 
during which the flux remains constant within statistical uncertainties. The time windows are enumerated 
at the top of the panel. The corresponding flux is shown by the blue marker below each number. The SED 
for each of the time windows is shown in Figure 6. 


flux hypothesis. The algorithm to determine the 
optimal partition is described by Jackson et al. 
(2005). The BB-binned light curve is shown in 
Figure 5. It is statistically compatible with the 
original light curve ( Xr/ndf — 257/232). This im- 
plies that flux variations within each BB cannot be 
distinguished with confidence from a locally con- 
stant flux. The shortest BBs are detected at the 
maximum of both sub-flares and have durations 
of ~9 hours. 

In order to measure the rate of flux increase at 
the rising edges of the sub-flares we parametrized 
them with an exponential function plus a constant 
background. The best-fit functions are shown in 
Figure 5. The time ranges over which the fits were 
performed were defined by the centers of the BBs 
before and at the maximum of each sub-flare. The 
resulting doubling time is 4.0 ±1.0 hours and 7.0 
±1.6 hours for the first and second sub- flare, re- 
spectively. As these values depend on the some- 
what arbitrarily chosen parametrization and fit 


ranges, we conservatively estimate that the dou- 
bling time scale in both sub-flares is < 8 hours. 

The PDS of the April 2011 flare is shown in 
Figure 4. It was obtained by computing the 
Fourier transform of the autocorrelation function 
using an algorithm for unevenly sampled data 
(Edelson & Krolik 1988). The PDS can be de- 
scribed by a power law of index « 1 . 1 and reaches 
the noise floor at a frequency of «0.6 cycles per 
day. The doubling time of the corresponding sinu- 
soidal component is s»10 hours, in agreement with 
the expectation from the measured doubling times 
of the flares. 

The pulsar flux remained unchanged during 
the flare, with an average flux above 100 MeV 
of Fp — (21.7 ± 1.1) x 10" ' cm -2 s _1 dur- 
ing the main part of the flare (MJD 55663.70- 
55671.02). The flux increase is phase-independent. 
This is illustrated in Figure 1. where the pha- 
sogram during the main flare period is shown. 
The peaks in the on-pulse interval remain at the 




Energy [MeV] 

•to 2 igf icr icf io 2 icf icf icf 


5 

N 

\ 

V. \ 

*****%, 

s 

T 

V, 1 \ 

\ 

v : 

V \ 

\ 

\ 

N 

\ \ 
\ 

1 \ V 

2 \ V" 

v \ 
3 . ■' 

V \ 

4 \ 

fry 

\ v 

' \ 

\ V 

\ \ 

\ 

\ 

\ N 

V 

' \ 

\ 

' \ 

\ \ 

\ 

V \ 

\ 

\ Vr 

\ \ 

• \ 

v \ 

V 

\ 

5 \ \ 

e \ V 

7 \ ' 

\ \ ..... 

8 \ 

c’ "h 

fy 

V. \ 

V \ 

\ ■ 

\ \ 

\ 

V \ 

'V 


9 \ ' 

10 '• V 

ii NA 

1 _ 



Fig. 6. — Spectral Energy Distribution evolution during the April 2011 Crab flare. The time windows are 
indicated in the bottom left corner of each panel and correspond to the ones indicated in Figure 5. The 
dotted line shows the SED of the flaring component, the dot-dashed line the constant background from the 
synchrotron nebula, and the dashed line is the sum of both components (see text). The average Crab nebular 
spectrum in the first 33 months of Femii observations is also shown in gray for comparison. 


same position. We also searched for periodici- 
ties other than the Crab pulsar with the time- 
differencing technique (Atwood et al. 2006), ap- 
plying the event-weighting technique described in 
Bickel et al. (2008). We scanned the frequency- 
range 0.1-256 Hz, allowing for a possible spindown 
up to twice the value of the Crab pulsar. No sig- 
nificant signal was found besides the pulsar, which 
was detected with a significance > 5.5<r. Finally, 
we searched for photon clumping on time scales 
shorter than the «10 min time binning by apply- 
ing a Bayesian Block analysis on the single photon 
arrival times, with no significant detection. 


4.1. Spectral evolution during the flare 

In order to measure the energy spectrum dur- 
ing the flare, and its evolution with time, the data 
must be averaged in time intervals long enough 
to ensure adequate photon statistics, but short 
enough to provide adequate temporal resolution. 
The 11 bins of approximately constant flux, de- 
rived from the BB analysis, provide a reasonable 
compromise between these two constraints. 

The SEDs for each of the time bins are shown 
in Figure 6, after subtracting the steady emis- 
sion from the pulsar and the inverse-Compton 
component of the nebula. It can be clearly seen 
that a new spectral component emerges from the 
synchrotron nebula during the flare, moving into 
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the Fermi energy range as the flare evolves. Its 
flux reaches a maximum between MJD 55666.997- 
55667.366 (frame 7); during this period the peak 
in the SED is clearly detected at E pea k = (375 ± 
26) MeV. 

It is difficult to parametrize the spectral shape 
of the flaring component due to the likely contam- 
ination by background flux from the synchrotron 
nebula not related to the flare. The determination 
of the latter is degenerate with the measurement 
of the flare component, as only the summed flux is 
measured. To break this degeneracy we proceeded 
under the following assumptions: 

1. The spectrum of the synchrotron nebula 
during the flare can be described by a power- 
law function and does not vary in time. 

2. The spectrum of the flaring component can 
be described by a power-law with an ex- 
ponential cutoff (equation 3 with k = 1). 
While the cutoff energy and normalization 
of the spectrum vary, the spectral index re- 
mains constant during the flare. 

We derived the spectral index of the flaring com- 
ponent and the spectrum of the background syn- 
chrotron component in a composite likelihood fit 
to all the time window's displayed in Figure 5, si- 
multaneously measuring the energy cutoff and flux 
normalization evolution of the flaring component 
in each of the time wfindows. For this we used 
the composite likelihood 2 part of the Fermi 
Science Tools. 

The best-fit values for the background syn- 
chrotron nebula during the flare period are Fs = 
(5.4 ± 5.2) 10“ ' cm -2 s _1 and 7 s = 3.9 ± 1.3. 
consistent with the average value measured during 
the first 33 months of observation. The spectral 
index of the flaring component is measured to be 
7 f = 1.27 ± 0.12. The best-fit values for Ep, c and 
the energy flux above 100 MeV, fp , are shown in 
Figure 7. 

This simple parametrization yields a good de- 
scription of the flare evolution, as can be seen in 
Figure 6. The data are therefore compatible with 
the interpretation that a new? spectral component 
of a fixed spectral shape evolves during the flare. 
As shown in Figure 7 the cutoff energy of the flar- 
ing component varies significantly, having a statis- 
tical probability of being constant of only 0.04%. 
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Fig. 7. — Total energy flux above 100 MeV, fp, as 
a function of cutoff energy E c for the flaring com- 
ponent of the April 2011 flare. The values corre- 
spond to the parametrizations shown in frames 1- 

10 of Figure 6 (the values obtained from frame Nr 

1 1 are not included, as no significant spectral cur- 
vature was detected in his time interval, allowing 
no robust determination of E c ). The number next 
to each marker denotes the corresponding frame. 
The axis on the right hand side indicates the ap- 
parent luminosity in units of the pulsar spin-down 
power of 5 x 10 38 ergs s -1 (Hester 2008). The 
numbered solid lines indicated the slope of the 
corresponding power-law dependency fp ~ Ep <c a . 
The dotted line marks the best fit function in this 
parametrization with a = 3.42 ± 0.86. 

However, while fp varies by more than an order 
of magnitude, the cutoff energy varies only by ap- 
proximately a factor of two. The relationship be- 
tween both quantities can be described approxi- 
mately by a power-law function fp ~ Ep , c a , with 
best-fit index of a = 3.42 ± 0.86. 

5. Discussion 

A year after first being reported, the gamma- 
ray flares from the Crab remain enigmatic. Where 
within the nebula does the emission come from? 
What produces the flux variations? How were the 
emitting particles accelerated? Ilow r are the flares 
related to the variability observed on yearly and 
monthly time scales? Although several ideas have 
been proposed, no certain answers can be given 
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today. The observations presented here give us 
the most precise look into the flare phenomenon to 
date, during the brightest outburst detected so far. 
We will proceed to discuss some of the implications 
and challenges posed by these observations. 

One striking property of the Crab nebula flares 
is their rapid flux variations, doubling within tj < 
8 hours at the rise of the 2011 April flare. Causal- 
ity arguments imply that the emission region is 
compact, with a length L < ctd « 2.8 x 10 -4 
pc. The emitted isotropic power at the peak of 
the flare of ss 4 x 10 36 erg s _1 corresponds to 
«1% of the total spin-down power of the pul- 
sar, the ultimate energy source of the nebula. 
It is difficult to explain how this energy is fo- 
cused into such a small emission volume. The 
focusing is generally easier to explain when the 
emission site is closer to the pulsar. The ab- 
sence of pulsation in the flare signal implies that 
the emission region is at least located outside the 
light cylinder of the pulsar. Another possibility 
to explain the flare brightness is that the emis- 
sion is highly anisotropic, as would be expected if 
the emission region moves relativistically toward 
us. While only mildly relativistic motion with ve- 
locities of «0.5c are observed inside the nebula 
(Scargle 1969; Hester et al. 2002; Melatos et al. 
2005), relativistic motion is expected in the pul- 
sar wind and in the downstream medium behind 
the wind termination (Camus et al. 2009). Rel- 
ativistic bulk motion is particularly expected at 
the “arch shock” of the wind termination, which 
has been proposed as the main site of gamma-ray 
emission (Komissarov & Lyutikov 2011). 

The flare emission is expected to result from 
synchrotron radiation by relativistic electrons and 
positrons (henceforth referred to together as elec- 
trons) (Abdo et al. 2011a). A new spectral com- 
ponent emerges in the SED during the flare. The 
hard photon spectrum (7 « 1.3) of the flaring 
component implies that most of the electron en- 
ergy is carried by the highest energy electrons. 
If the electron particle density 11 (e) per energy 
at an energy e is characterized by a power law 
n(e) ~ £~ p , the spectral index is p < 2j — 2 « 0.6. 
in a random magnetic field in which the electrons 
are isotropically distributed. The energy per log- 
arithmic energy interval, which is proportional to 
e x 11 (e), is therefore rising with energy. 

It is an interesting question how such a hard 


electron spectrum is produced. Standard diffu- 
sive shock acceleration typically results in spec- 
tra with p > 1 (Gallant et al. 1992; Kirk et al. 
2000). Even though it has been shown that harder 
spectra can be produced in certain field configura- 
tions with low-level turbulence (Kirk & Heavens 
1989; Summerlin & Baring 2011), these condi- 
tions are not expected at termination shocks of 
pulsar winds. Additionally, shock acceleration ap- 
pears to be inefficient at highly oblique shocks 
that are representative of the pulsar wind ter- 
mination discontinuity (Ellison & Double 2004; 
Summerlin & Baring 2011; Sironi & Spitkovsky 
2011b). One alternative is that magnetic re- 
connection in the striped pulsar wind might ac- 
celerate particles (Lyubarsky 2003; Kirk 2004; 
Yuan et al. 2011; Bednarek et al. 2011). How- 
ever, simulations show that reconnection behind 
the pulsar termination probably does not provide 
the required electron energies to produce gamma- 
ray emission (Sironi & Spitkovsky 2011a). An- 
other interesting possibility is that acceleration is 
occurring directly in the electric field induced by 
the pulsar, as discussed by Abdo et al. (2011a). 

The observation of a peak synchrotron energy 
of «380 MeV is among the highest yet seen from 
astrophysical source today. The observation is 
surprising as particle acceleration in the presence 
of synchrotron cooling is expected to limit syn- 
chrotron emission to photon energies below «150 
MeV (Guilbert et al. 1983; de Jager et al. 1996; 
Komissarov & Lyutikov 2011). Two solutions to 
this problem have been proposed recently in this 
context: 

1. The electric field at the acceleration site 
is larger in magnitude than the magnetic 
field. This is generally an unstable state 
in plasma, as charges will short out the 
electric field; however, temporarily such 
a configuration is expected, e.g. in mag- 
netic reconnection events. The gamma- 
rav emission might occur afterward, when 
the accelerated electrons enter a region of 
enhanced magnetic fields (Uzdensky et al. 
2011; Cerutti & Uzdensky 2011), More de- 
tailed studies are required to assess whether 
such a scenario is plausible in the nebula 
environment and can be sustained for the 
duration of the flares. 
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2. The gamma rays are emitted in a region 
of bulk relativistic motion, and are there- 
fore Doppler boosted toward the observer 
(Komissarov & Lyutikov 2011). For a flow 
moving directly toward us a Lorentz factor 
> 2 is sufficient to accommodate the ob- 
served peak energy. 

Variations in the Doppler boosting can nat- 
urally account for the observed flux variation 
(Lyutikov et al. 2011). The observed spectral 
evolution is compatible with such an interpreta- 
tion: the energy flux of the emission varies ap- 
proximately as a power of a = 3.42 ± 0.86 with 
the cutoff energy; a correlation with a ss 3 is in- 
deed expected for variations produced by changes 
in relativistic beaming (Lind & Blandford 1985). 

The flare brightness, the high frequency of the 
observed peak of the gamma-ray emission, and the 
spectral evolution during the flare all suggest the 
presence of relativistic beaming. We therefore con- 
clude that, independent of the location of the emis- 
sion region and the physical processes responsible 
for the flares, the emission region is moving rela- 
tivistically toward us, and changes in its motion 
are likely the predominant mechanism responsi- 
ble for the observed flux variations. Such a kine- 
matic explanation does however not address the 
issue of how a moving source can be created dy- 
namically and sustained radiatively in the face of 
strong losses. The Crab Nebula still has much 
more to teach us. 
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